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1. Introduction 

The finite-size method, originally proposed by Luscher [JTJ , allows us to calculate the scattering 
phase of two particles in infinite volume from the energy eigenvalues of a two-particle system 
enclosed in a finite spatial box. For a fixed finite volume, the scattering phase is determined at a 
discrete set of energies. To better capture the energy dependence of the scattering phase, we would 
like to increase the number of energy eigenvalues accessible to lattice calculations. Particularly, in 
a scattering channel where a resonance appears, we want to calculate the scattering phase at many 
energies in the resonance region in order to determine the resonance parameters such as the mass 
and width. 

One method to address this issue is to use different lattice sizes to obtain the scattering phase 
at more energy values. However, this requires additional simulations that are increasingly more 
demanding at large volumes. An alternative way is to use the moving frame (MF) technique, 
introduced by Rummukainen and Gottlieb [0], which generalizes Luscher 's original method from 
the center-of-mass frame (CMF) to a MF. The key point is that the energy spectrum calculated in a 
MF is different from the one obtained in the CMF. Thus combining the CMF and a MF allows us 
to compute the scattering phase at more energies using the same lattice size. This will increase the 
accuracy of the calculation of the desired resonance parameters for less computational cost. 

As an application, we have performed a calculation of the P-wave pion-pion scattering phase 
in the p -meson decay channel using three Lorentz frames: the CMF, the original MF with total mo- 
mentum P = (27T/L)e 3 (MF1) and a second MF with P = (2n/L)(e\ +e 2 ) (MF2). This approach 
allowed us to determine the scattering phase as a function of the energy covering the resonance 
region [Q]. MF2 is a new MF for which the finite-size formulae were not yet available in the liter- 
ature. In this work, we present a derivation of the finite-size formulae for this MF. Our derivation 
closely follows the work in Refs. [|l], [2p. 

2. Two-particle wave function 

Consider a two-particle system with total momentum P / 0. The state of such a system can be 
described by a wave function <j>(x), where x is the relative position between the two particles. To 
establish the formula for the scattering phase, which is naturally defined in the CMF, we need to 
transform the scattering system from the MF to the CMF. This is accomplished by the appropriate 
Lorentz transformation, parameterized by y, under which the wave function transforms from <j) to 

0CM- 

We assume that the two-particle interaction vanishes in the region where |x| > R, called the 
exterior region. Then in the exterior region (f>cM ( x ) satisfies the Helmholtz equation 

(V 2 + p 2 )0cm(x) = O, for |x| >7? , (2.1) 

where p is given by the energy-momentum relation p 2 = (Ecm/2) 2 — m 2 with Eqm the center-of- 
mass energy and m the single-particle mass. If we consider such a two-particle system enclosed 
in a box with finite size L > 2R, then the total momentum P is discretized as P = (27r/L)d with 
dGZ 3 and tpcM (x) satisfies the d-periodic boundary condition 

4>cm(x) = (-l) d >aw(x + ynL) , for all neZ 3 . (2.2) 
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The details of the derivations of Eqs. (2.1) and (2.2), along with any unexplained notation such as 
yn, can be found in Ref. 

A function that satisfies the Helmholtz equation (2T) and obeys the d-periodic boundary con- 
dition (2.2) is called a singular d-periodic solution of the Helmholtz equation. A simple example 
of such a solution is the Green's function 

keP A K P 

where the summation of the momentum runs over 

P A =Ue 



jkeM 3 k= ^7 _1 (m + id) , for m e I? J . 



More singular d-periodic solutions can be generated from the Green's function by introducing the 
harmonic polynomials 

ty m (x)=r l Y bn (±), r=\x\, (2.3) 

and defining 

Gt n (x,P 2 ) = m m (V)G d (x,p 2 ). 

In Eq. ( |2.3| ) the notation x represents the solid angle parameters (d, (p) of x in spherical coordinates. 
It was proved in Ref. [jl]] that the Gf m (x,p 2 ) are complete and linearly independent. Therefore 
0cm(x) can be expanded in terms of Gf m (x,p 2 ) as 

0cm(x) = 2> OT Gf m (x,/7 2 ) , for R < \x\ < L/2 . 

l,m 

As given in Ref. Gf m (x,p 2 ) can be expanded in terms of the spherical harmonics F/ m (x) 
and spherical Bessel functions ji (pr) and ni (pr) through 



Gf m (x,p 2 ) = ( -^fp !+l I Y lm {x)m{pr)+ £ ^n/Wl'neWhipr) \ ■ (2-4) 



Due to symmetry considerations, some ^f m [i m i(p) vanish. For the MF2 (d = ei +e2), we list all 
the non-zero values of ^f m Vm i (p) for 1,1' = 0,1 in Table [jj There we expand ^f m Vn j (p) in terms 
of the modified zeta function, which is defined by 

, _ P L 

in. 



(\n\ 2 -q 2 Y' * 2?r 

This zeta function is formally divergent and needs to be analytically continued. Ref. |§] gives a 
numerically calculable expression of ^ n (\\q 2 ), which is not, however, valid for the MF2. Here 
we give a more general expression that can be used in all the different MFs 

1 „ wvn ir . „ !^7U| 2 



Jo uezwo ? f 1 

+ 7 fdtiff* - i)J^8 I0 d m{ A 3/2 - yn8 l0 8, 
Jo \/4k t 



£ T §^exp(-(|n| 2 -^ 2 )). 

^neP d q 
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ll 


-^^(l;^ 2 ) 


ll 
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Table 1: ^f m expanded in terms of ^J^(l;^ 2 ) for /,/' = 0, 1. The notation m = —to is used. 



3. Deformed symmetry 

We introduce the group G as the set of all lattice rotations and reflections R that leave the set 
of Pd invariant, 

G = {R\Rk£P d , VkeP d }. (3.1) 

Let J$f be the space of all d-periodic wave functions y(x). For any i//"(x) G and fieG, the 
transformed wave function Ry(x) = y(R~ l x) satisfies 

¥ (R- l x) = £ Vlm m m (R- l V)G d (x,p 2 ) = £ v,nA{R)%n'{V)G A {x,p 2 ) G ^ . (3.2) 

l.m l,m,m' 

This shows that J4? (or equivalently the two-particle system) is closed under the group G. i (R) 



used in Eq. ( |3.2| ) is the standard finite-dimensional rotation matrix. It originates from 



m m (x)=m,n(R- i x)= £ D2Lr(^)9W(x) . 

m'=-l 

If we consider a vector space Jtf in which the harmonic polynomials &i, n (x) form an orthonormal 
basis, then D^ m ,(R) is simply an irreducible representation (irrep) of the rotational group 0(3), 
which describes the group elements R G 0(3) in terms of linear transformations acting on the 
vector spaces Jfi. If we restrict R to R G G, then , (R) is reducible and the vector space M\ can 
be further decomposed to sub-spaces that are invariant under G. 

In the case of d = (7 = 1), G is given by the full cubic group Oh- For the MF cases, the 
constraint in Eq. ( |3.lD excludes some lattice rotations and G is reduced to a subgroup of Oh- One 
can prove that Eq. ( [Q| ) is equivalent to 

G = {R g O h I Rd = d or Rd = -d } , 

indicating that for d 7^ G is the parity doubled little group of P = (2n/L) d. This is the group of 
rotations that leave the specific direction d unchanged combined with reflections of d. 

In the MF2, G is given by the orthorhombic group D211, which has 8 one-dimensional irreps: 
A ± , Bf, B2 and B*. 1 The index ± comes from parity, which fixes the transformation behavior of 
0cm (x) under reflections x — > —x. 

'With the change of coordinates = +X2), x 2 = ^i* 1 ~ x z) an d x '\ = x 3> tne notation of the irreps used 

here coincides with the ones used in Ref. ftm, chapter XII. 
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In general, we denote the irreps of G as T. As we mentioned, the space Jtf, which is invariant 
under 0(3), can be decomposed into the irreducible sub-spaces J#r of D2/,. For brevity, we rep- 
resent these decompositions by only the indices of / and T. For / = 0, 1,2 the decompositions are 
given by 



/ = decomp A + 
1=1 decomp B\ 
1 = 2 decomp A + 



3 > 



The corresponding basis polynomials of J^r can be written as \T,a;l,n), where a runs from 1 to 
dr, the dimension of T, and n runs from 1 to N(r,l), the total number of occurrences of T in the 
decomposition of the space J^j. The relation between \T,a;l,n) and ^/ m (x) is written as 

/ AT(r,Z) d r 

\F,a;l,n) = £ <£ iew »3£*(x) , ^ m (x)=£ £ £ cft,a,n;ml r 5 a;l,n) , 



m=—l 



T n=l a=\ 



where the coefficients c' r a n . m form a (2/ + 1 ) x (2/ + 1 ) matrix C\ satisfying Cj C\ = 1 . We list 
|r, GC;l,n) in terms of 2^ m (x) for / = 0, 1 in Table ||. 



r,a;/,n) 


in terms of Wi m (x) 


in terms of polynomials 


|A+,1;0,1) 


%o(x) 


1 


l*T,i;i,i> 






|B 2 ,1;1,1) 


-^fu(x) + H^(x) 




1,1) 


^io(x) 





Table 2: |r, a;l,n) expanded in terms of £%„(x) for / = 0, 1. The notation m = —to is used. 



The matrix (R) can be diagonalized through 

/ D (r0 
QD^C] = £> (r2) 

V 

where D( r \R) is the matrix representation of Dih for the irrep T. It gives the action of the group 
element R G G in terms of linear transformations on the vector space J£f as follows 

R\T,a;l,n) =D ( p p {R)\T,p;l,n) . 



From inspecting Eqs. ( p.4| ) and (p.2|), it can be seen that the matrix ^ , transforms as 



According to Schur's lemma, we then have 



C, 



T.a,v.m- " Im.l'iii 



C r\a'n':m' = ^r,V^a,a'Mf Vnl {T) . 



(3.3) 
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Using the projection operator 



G ReG ReG 



we can project (pcM ( x ) to the T representation through 0^(x) = Pr.a<PcM( x )- Following from the 
orthogonality theorem 



one can prove that 



N(T,l) 



^Cm( X ) ~ H ^1 V lm c r,a,n;m c T,a,n-jn'Gl m '(x,p 2 



l,m,m' n=l 



J>„ \ r- l \T,a^n)n l (j>r) + Y, M fnM T ) r \r,a;l',^)j P (pr) \ (3.4) 

l,n \ l',n' J 



where the coefficients V/„ are defined as 



Vln 



An 



m=-l 



4. Finite-size formulae 



In Eq. (|3~4|), we have expanded 0^(x) in terms of the singular d-periodic solutions of the 



Helmholtz equation. In the exterior region, it can also be expanded in terms of the spherical har- 
monics Y[ m (x) and spherical Bessel functions j) (pr) and ni (pr) as 



<t>CM ( x ) = L b lm Y lm (±) (<% (p)ji (pr) + (p)n t (pr) ) . 



(4.1) 



l,m 



It was proved in Ref. [Ill] that there exists a unique eigenfunction of the Hamiltonian in the infinite 



volume which coincides with 0™ (x) in the exterior region. As a result, the coefficients CCi(p) and 



fii(p) in Eq. (4.1) can be related to l-th wave scattering phase 8i(p) through 



w 2/5,(p) = a '(p) + i P'(p) 
oci(p)-iPi(p) 



Eqs. (41) and (^4) together determine the coefficients b\ m that make the two expansions of 0^(x) 
coincide. Writing bi m as b\ m = Y, n bi„c l r a trm , Eq. ( jO] ) can then be written as 



0cm( x ) =Y) ) ^ r '\r,a;l,n) (a,(p)ji(pr) + p,(p)n,(pr)) 

l,n 



We therefore get the relation between V/„ and bi n as 



bi n pi(p) = Vi,i , 

~bi n a,(p) = 1 £vi'n'Mf n , ln (r) . 



I'n' 
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This homogenous system has a non-trivial solution of bi„ only when the determinant of the coeffi- 
cient matrix equals zero, which results in 

det [§ -M d (r)] = , &„, /v = fydnj tan" 1 S,(p) . (4.2) 



If we consider only the lowest angular momentum contribution to Eq. ( |4.2[ ), then the finite-size 
formulae are given by 

tan- 1 5 (/>)=M % 1 (A+), 

tan" 1 ft (p) =M? M1 (r) , for r = B7, B 2 , . 
Using Tables [l], ^ and Eq. (p.3|), we can describe Mf n Vn ,{T) in terms of ^(Uq 2 ) 

M A om {A + ) = (Y^ ,2 qr x ^ 

<n(fiD = (Yn 3/2 q)-\^ ~ ^^ + i^=^(2?£ - iF 2 d 2 )) 

Mfi.n^) = (y* 3/ \rWo+^^o) • 

Thus we obtain the finite-size formulae in the MF2 for the irreps T = A + , B^, B^ and B^. 

5. Conclusion 

The finite-size formulae for a new MF with total momentum P = (2n/L)(ei +^2) are derived 
in this work. These formulae can be used to calculate the S-wave and P-wave scattering phases. 
Using similar procedures to those in this work, one can easily generalize the derivation to other 
MFs and different irreps. 
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